Pseudo-Quantum Criticality in Electron Liquids Exhibited in Expanded Alkali Metals 



Hideaki Maebashi and Yasutami Takada 
Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, 

(Dated: February 1, 2008) 



Japan 



o 
o 

(N 
G 



5h 



i 

C 

O 

o 



o 
o 

o 
o 



X 



With paying special attention to the divergence in the compressibility k, we study the Coulombic 
screening in alkali metals to find singular long-range fluctuations in the electronic polarization 
originating from this divergence. As a consequence of this singularity, we predict the decrease of 
the equilibration distance between ions against the increase of r 3 the Wigner-Seitz radius of valence 
electrons, provided that the condition of 2r c < r s < 4r c is satisfied with r c the ion-core radius. This 
prediction is in good quantitative agreement with the recent experiment on liquid Rb. 
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According to the recent high-resolution measurement 
of x-ray diffraction for expanded liquid Rb [l[ , an average 
distance between nearest- neighbor ions i?NNi a $ deter- 
mined by the first-peak position in the radial distribution 
function g(R), decreases with increasing Ro the Wigner- 
Seitz (WS) radius of ions. This is in sharp contrast with 
the normal situation in compressed liquid alkali metals in 
which i?NN increases with Rq 0, H, 0, 1H . In solid Rb the 
valence-electron density n is specified by r s — 5.2 with r s 
being the WS radius of electrons (in atomic units which 
we take throughout this Letter), very close to r s — 5.25 
at which the compressibility k of the 3D electron gas di- 
verges |H . Thus n in expanded liquid Rb is in the range 
of r s > 5.25 where k is negative |7|. In this Letter, we 
shall explain this unexpected contraction of -Rnn by fo- 
cusing on the divergence, rather than the negativeness, 
of k as a key issue. 

At arbitrary n, the compressibility sum rule re- 
lates k to the polarization function H R (q,u>) (with 
the superscript R denoting the retarded form) through 
lim 9 ^o n fl (g, 0) = n 2 n. Hence, near the critical density 
of r s — 5.25, H R (q,uj) becomes singular in the limits of 
q/pF —> and uj/vpq — > with pp and vf being, re- 
spectively, the Fermi wave number and the free-electron 
Fermi velocity. More specifically, we find that H R (q,uj) 
is well characterized by such a singular form as 



n«( g , w ) 



+ q 2 -icj/Tq), (1) 



where £ 2 = T = (2/tt)« p £ f , and £ F 2 = Up 2 

with Kp the free-electron compressibility. Note that near 
the quantum critical point (QCP) of dynamical expo- 
nent 3 [H, the response function takes exactly the same 
form as in Eq. (fTJ) . In the electron gas, however, there 
is no true phase transition (and in fact "the correlation 
length" £ turns into a pure imaginary value for negative 
k), owing to the fact that H R (q, u>) in this case represents 
response to a total of external and induced charges, mak- 
ing it acausal . One of the interesting consequences of 
Eq. ([T]) near this "pseudo-quantum critical point" is the 
softening of excitonic collective modes, as featured by the 
dispersion relation of u> = (2/w)\ Kf/ k + q 2 £ F \vFq [10|]- 

The classical or Coulombic screening of ions near this 
pseudo-QCP poses a distinct problem from the quan- 
tum or Kondo screening near the true QCP and 



deserves special attention, basically because it is useful 
in clarifying the important physics of the coupling be- 
tween the q = singularity in the Coulomb interaction, 
v(q) = 47r/(7 2 , and the singular fluctuations in U R (q, 0). 
Another important aspect of the problem is the fact that 
the ions in alkali metals are not point charges but com- 
plexes of a nucleus and core electrons surrounding it, the 
feature of which can be well captured by a suitable pseu- 
dopotential V ps (q). This complex structure combined 
with the pseudo-QCP makes the observable phenomena 
intriguing, in which an important role is played by the 
presence of a characteristic wave number qo (proportional 
to the inverse of the core-electron radius r c ) at which 
Vp S (q) vanishes due to the orthogonality between valence- 
and core-electron wave functions [l2j . 

In this Letter, we report on our finding that the 
pseudo-quantum criticality is intimately connected with 
the appearence of a local minimum in the electrostatic 
potential <p(q) at q = for some range of r s (see Fig.[lJ. 
This minimum is shown to be brought about by (i) the 
Thomas-Fermi-type screening with the singular polariza- 
tion and (ii) the subsequent rearrangement of the charge 
distribution in response to the resulting screening hole. 
The first process induces the contraction of the equilibra- 
tion distance between ions, while the second acts in the 
opposite way. In expanded alkali metals, this second pro- 
cess is suppressed by the presence of go in V ps (q), leading 
to the experimental results [l[ with which our theoreti- 
cal results of g(R) obtained by systematic Monte-Carlo 
simulations for Rb ions are in good accord. 

In terms of H R (q,0), we can write 4>{q) induced by a 
point charge in the 3D electron gas as 

<f>(q)=v(q)/[l + v(q)Il R (q,0)}. (2) 

Using the exchange-correlation kernel f^ c (q,uj), we 
can give a formally exact experession for H R (q,ui) as 
U R (q,cj) = \n? R (q,u)- 1 +fxc(q,u)]-\ where U 0R (q,cj) 
is the free-electron polarization function. Since an ac- 
curate parameterization of the diffusion Monte Carlo 
(DMC) data for f xc (<l,0) is available [13j, the virtually 
exact 4>(q) is known. (See the solid curves in Fig. [lja).) 
The DMC data reveal that f xc (q,0) is flat near q = 0, 
e.g., / xc (q,0)-lim^oixcM) ~ 0(q w ) for r s = 5 0, 
allowing us to derive the result in Eq. ([T]) for q/pF 1 



2 



sfs=1 



(^(q)/v(q s ) - 
* (s) (q)/v(q s ): 



r s =3 



r s =1 1 



(b) 



(p F / q )n W (q/p F ) 



r s =5.25 
T\\ (5=0) 



x 



i \ 



r s =3.17 



w 

\r s =8.43 

¥ s r 1) 



0.8 



0.6 



0.4 



0.2 



0.0 



0.5 1.0 

q/q s 



.50 



1 2 3 
qr s /7i 



o.o 



FIG. 1: (a) Electrostatic potential <f>(q) scaled by v(q s ) with 
Qs = \J 9tf /£f as a function of q/q B with changing r s from 1 to 
11 by 1 (the solid curves from top to bottom), together with 
the corresponding screening potential 4>^ s ' (q) in Eq. Q by the 
broken curves, (b) Rearrangement of the charge distribution 
due to the screening hole n' r '(g/pF) divided by q/pF, which 
is properly scaled in Eq. (J6j , as a function of qr s /ir. The 
"distance" from the pseudo-QCP is defined as S — </>(0)/v(q s ). 



and Lo/vFq <C 1 under the assumption that the limits of 
q —>■ and u> — ► are interchangeable for / xc (g, w). 

Although H 0R (q, 0) itself is specified only by pp = ar s 
with a = (4/97T) 1 / 3 , a unique length scale for this func- 
tion, consideration of the coupling to the Coulomb inter- 
action evokes another length scale, g^p = {■Kar s /A) 1 l 2 , 
as characterized by v(q)U 0R (q 7 0) ~ q\yjq 2 for q — > 0. 
For the full polarization function n fl (g, 0), Eq. ((T|) indi- 
cates |£| as a characteristic length scale which becomes 
much longer than £p ~ P-p 1 for k/kf 1. Eventually at 
the pseudo-QCP, concomitantl y with t he divergence of 
|£|, a new length scale, q^ 1 = y/^p/ qrF , emerges for the 
coupling to the Coulomb interaction, as characterized by 
v(q)n R (q,0)~qt/q 4 forq->0. 

Keeping those length scales in mind, we represent (j>(q) 
in an expansion form around q = for arbitrary n as 



cb(q)=v(q s )[S+(l-5 2 )(q/q s ) 2 + 0(q i )] 



(3) 



Here we have introduced 6 [= (kf/k)/(<Ztf£f)] as a key 
quantity describing the "distance" from the pseudo-QCP. 
In Eq. ([31), the coefficient of the q 2 term is positive for 
\5\ < 1 corresponding to 3.17 < r s < 8.43, which clarifies 
the reason why 0(g) in Fig.rjja) exhibits a local minimum 
at q = only in this range of r s . We shall call this 
range the pseudo-quantum critical region, in which alkali 
metals, Na, K, Rb, and Cs, are included. 

In order to obtain a deeper insight into the Coulombic 
screening in the presence of the singular contribution in 
H R (q,0), let us write <j>(q) as 

0(g) =«(«)[ 1 - n ind (g) ] - ^ (q) + 0« (g) , (4) 

where n ind (q) [= v(q)U R (q, 0)/(l + v(q)U R (q, 0))] is the 
electron density distribution induced by the point charge 



and 0( s ) (q) is defined by substituting the singular polar- 
ization, Eq. |1]), into Eq. © as (/)^(q) = v(q)[l-n^ d (q)] 
with n£ d (q), given by 



n^(q) = [l + S(q/q s f + (q/q s r 



(5) 



It is easily verified that when r s — * 0, Eq. ([5]) is re- 
duced to the standard Thomas-Fermi screening charge, 
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(?) — * <?tf/(9tf + Q 2 )- Thus we may regard </>^(g) 
as the Thomas-Fermi-type screening potential with the 
singular polarization. In Fig. QJa), <P (q) 1S plotted by 
the broken curves with changing r s from 1 to 11 by 1. 

The residual term in Eq. (0]), ^^(q), is the electro- 
static potential due to the rearrangement of the induced 
charge from n^ d (q) to n- md (q), amounting to 4>^(q) — 
v(q)n( T \q) with n^\q) = n^ d (q) — n in ^(q) . Since the sin- 
gular long-range part n^Jq) is extracted from ni n( j(g), 
we speculate that (q) is characterized only by the 
length scale pp 1 , implying that n^(q) should be ex- 
pressed in terms of some universal function of q/pF- This 
speculation is checked by casting rv- T '(q) into the form as 

n (r) (?) = (2 + Sr'iQTF/qs) 2 nW(ff/pp) • (6) 

Here we have determined the prefactor by respecting the 
result of n^ d {q s ) = (2 + S)- 1 . In Fig. C^b), we plot the 
obtained fS r \x)/x with x = q/pF for (5 = and ±1, from 
which we see that h^\x)/x, having a peak structure for 
k/t s < q < 27r/r s , may well be regarded as universal in 
the entire pseudo-quantum critical region. 

In calculating Va(R) the effective interaction between 
ions, we resort to second-order perturbation with respect 
to V ps (q) [l4|]- For R outside the core region, we obtain 
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Pps(q) 2 4'(q)qsm(qR)dq, (7) 



where p ps (g) is the effective charge of the ion, given by 
Pp S (q) = (g 2 /47r)Vp S (g). In general, the charge neutrality 
of the whole system imposes p ps (0) = 1. For r c > 0, 
Pps(q) vanishes at q = go- I n the limit of r c — > (or 
go — > oo), however, p ps {q) — 1 for arbitrary g, which re- 
duces Vii(i?) into (j>(R), the inverse Fourier transform of 
4>(q). In the following specific calculations, for simplic- 
ity, we employ the Ashcroft empty-core pseudopotential, 
for which p ps (q) = cos qr c , but our conclusions do not 
depend on the detailed shape of Vp S (g) very much. 

In accordance with Eq. we shall devide Vn(R) in 
Eq. into two components as 



Vn(R) = v£ s) (R) + vy>(R) 



where V» B \R) and Vj; l; (i?) are, respectively, defined by 
using (/)^(q) and <jft\q) instead of </>(q) in Eq. ([7]). Since 
n|^ d (g) is given in a simple analytic function, we can an- 
alyze (R) by pursuing the poles of <fi^ (q) or equiva- 
lently those of ^S(g) in the complex g space. According 
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FIG. 2: Effective interaction between ions, Vu(R), with using 
the empty-core pseudopotential with r c = 2.0, for r s = 5, 6, 
and 7 (the solid curves). The broken curves represent V^IR) 
in Eq. <(9j . The inset shows the trajectory of the screening pole 
q a e lBt> with increasing r a in the complex q plane. 



to Eq. ([5]), all the poles are on the imaginary axis for 
5 > 2, while they deviate from the axis for 6 < 2 (see the 
inset of Fig. [2]) , implying that the "screening length" in 
the Thomas-Fermi-type screening potential changes into 
a complex number, so that the resulting potential shows 
a damped oscillatory behavior. As 6 decreases, the oscil- 
lation amplitude increases and eventually for S < 1, the 
short-range attractive part develops in VJj (.R). Note 
that this oscillation related to the pseudo-quantum criti- 
cality is distinct from the Friedel oscillation, because the 
latter does not originate from those poles but the branch- 
cut singularity in <p(q). 

For 1 5 1 < 2, the poles in the upper- half complex q plane 
reside at q s e Ws and -q s e- w * with 6 S = (1/2) cos _1 (-<5/2). 
Writing ,Op S (g s e l9s ) as Z ps e ie>ps , we obtain vS s \r) as 



i/OO/di - 72 sm(2ip - q s R cost 
V * {R )- Z v* i?sin20 s 



-q s R sin t 



ps 



(9) 



with ip = 9 S — Op S + 7r. This analytic result with r c = 2.0 
is plotted in Fig. [2] by the broken curves for r s = 5, 6, 
and 7, in comparison with the corresponding numerical 
result for Vn(R) (the solid curves). The global minimum 
in Vn(R) determines the equilibration distance between 
ions, R cq , which is found to decrease with increasing r s . 

(s) 

The same is true for Rcq the equilibration distance deter- 
mined through V~ s \r), suggesting that the main feature 

(s) 

of i? oq will be studied by examining R cq . An approxi- 
mate result for R co ] is obtained from Eq. (J9j) as 

Rg> a (2(^ + 7r/2)/(g s cos(9 s ) = (Ap + 7r)/q s V2~S . (10) 

A rather extensive examination of various factors in 
Eq. © reveals that R cq is mainly controlled by tp. As 
shown in the inset of Fig. [21 the phase 9 S decreases as the 



pole q s e l6s evolves along the trajectory with increasing 
r s . On the other hand, the phase ps increases, because 
it is approximately given by the phase of q — qo itself, as 
can be seen by expanding p ps (q) around q = qo- There- 
fore, in a combined manner, both changes contribute to 
the decrease of the phase <p, leading eventually to the 

(s) 

contraction of R cq with the increase of r s . Incidentally, 
for very large qo as in the case of a point charge, ps 
hardly changes. Thus in such a case, only the phase 9 S 
contributes to the decrease of R C q ■ 

The presence of VJj (i?) in Eq. ([5]) manifests itself in 
the difference between R cq and R oq . This residual in- 
teraction is determined by n^ T \q) in Eq. ©. Since it 
is scaled by Pp 1 , the resulting interaction is necessarily 
scaled by the same length scale or ar s , indicating that 
V^ T \R) contributes to increasing R cq with r s . In this 
context, it is important to note that n^(q) has a sharp 
peak structure as shown in Fig[T] (b). Thus, if p ps {q) is 
small in the peak region of n^ r ' (q) as realized by the con- 
dition of n/r s < qo <2n/r s , the contribution of V» T \R) 

(s) 

is strongly suppressed, making R cq approach R cq . This 
suppression never occurs in the point-charge case. 

In Fig. [3l we show the overall behavior of R oq with 
the change of r s /2r c , which is obtained numerically with 
using the empty-core pseudopotential. The broken lines 
correspond to 5 = 1 (r s = 3.17), S — (r s — 5.25), and 
5 = —1 (r s = 8.43), showing the pseudo-quantum critical 
region. The contraction of R cq is clearly seen for 2r c < 
r s < 4r c in accordance with 7r/r s < qo < 27r/r s , where 
qo = 7r/2r c . Although the global minimum of Va(R) 
is obtained at its first minimum for 6 < 1, the second 
minimum becomes lower than the first one for S > 1, 
making R cq jump, which is not shown in Fig. [3l For 
S < — 1, the situation is completely different because of 
the shift of the global minimum of 4>{q) to a finite q (see 
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FIG. 3: Equilibration distance between ions, i? cq , for the 
empty-core pseudopotential with changing r c from 1.0 to 3.0 
by 0.1 versus the ratio of the WS radius of the electrons r s to 
the core diameter 2r c . 
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FIG. 4: Comparison between our theoretical results of g(R) 
obtained by the canonical Monte-Carlo simulations (the open 
circles) and the experimental results in Ref. [1] (the solid 
curves). The broken line denotes -Rnn of g(R) at 373 K. 

Fig. [1]) . This may be related to the instability toward the 
Wigner transition. 

The contraction of R cq may be explained from a geo- 
metrical point of view. Inserting two ions into the elec- 
tron liquid excludes the valence electrons from the core 
regions owing to the orthogonality between valence- and 
core-electron wave functions. For the case of 2r c < r s < 



4r c , the rearrangement of the induced charge distribution 
due to the screening hole is restricted, because the WS 
radius r s , which is the typical length of the electron liq- 
uid, is comparable with the diameter of the excluded vol- 
ume. Then, Va(R) is dominated by the complex screen- 
ing V» b \r) in the pseudo-quantum critical region. On 
the other hand, when r s > 4r c , the valence electrons rec- 
ognize the ion core virtually as a point charge, so that 
the problem is reduced to R cq between two point charges 
which is an increasing function of 7' s . In Fig. [3l we have 
also plotted R eq for solid alkali metals, Na, K, Rb, and 
Cs. Since r s /2r c ~ 1 in these metals, the contraction of 
i? eq is generally expected in the expanded alkali metals. 

Finally, in Fig. |4j we present the canonical Monte- 
Carlo results of the radial distribution function g(R) for 
500 classical particles interacting with each other through 
Va(R) in Eq. ||7J) at finite temperatures, together with the 
experimental results in Ref. For liquid Rb, we have 
used the empty-core pseudopotential with r c = 2.4 [l5| . 
The agreement between theory and experiment is good 
not only for the first-peak position but also for the overall 
change of g(R) with increasing r s and temperature. 

In summary, we have shown that the pseudo-quantum 
criticality in electron liquids transforms the standard 
Thomas- Fermi screening into the complex screening char- 
acterized by the development of a short-range attractive 
part in the screened Coulomb interaction. For 2r c < 
r s < 4r c , the core- valence orthogonality promotes this 
complex-screening effect by suppressing the rearrange- 
ment of the induced charge due to the screening hole, 
leading to the contraction of the equilibration distance 
between ions with increasing r s in good agreement with 
experimemt on expanded liquid Rb. 
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